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Abstract 

Attempts to apply effective field theory (EFT) methods to nonrelativistic 
nucleon-nucleon (NN) scattering have raised questions about the nature and 
limitations of an EFT expansion when used nonperturbatively. We discuss 
the characteristics of a meaningful EFT analysis and compare them with tra- 
ditional approaches to NN scattering. A key feature of an EFT treatment is a 
systematic expansion in powers of momentum, which we demonstrate using an 
error analysis introduced by Lepage. A clear graphical determination of the 
radius of convergence for the momentum expansion is also obtained. We use 
these techniques to compare cutoff regularization, two forms of dimensional 
regularization, and the dibaryon approach, using a simple model for illus- 
tration. The naturalness of the parameters and predictions for bound-state 
energies are also shown. 



I. INTRODUCTION 



The study of nuclear forces and nucleon-nucleon (NN) scattering has a long history [|TJ . 
Numerous models have evolved that provide superb reproductions of data over large energy 
ranges with only a few parameters per scattering channel. Nevertheless, the models do not 
provide reliable error estimates, and their connection to the symmetries and dynamics of the 
underlying theory of the strong interaction, quantum chromodynamics (QCD), is obscure. 

Recently, many papers have addressed the application of effective field theory (EFT) 
methods to NN scattering The hope is that such an analysis will yield a system- 

atic, controlled expansion and bring new understanding to how various processes involving 
nucleons relate to each other This requires a nonperturbative treatment of the EFT to 
properly account for the interaction between two or more heavy nucleons ||. 

There are currently many disagreements in the literature about the nature and limitations 
of an EFT expansion in this case |5|-|5,n]. Regularization is required to handle divergences 
that arise, but the results are said to depend on the regularization scheme used (see Table [[[) 
and the size of the scattering length involved. More generally, it is claimed that the behavior 
and predictive power expected from a true effective field theory is not exhibited by every 
regularization method when applied nonperturbatively |5|JJ. However, others claim that 
results are independent of the regularization method when proper power counting is applied 
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Our goal here is to clarify the important features of a nonperturbative EFT. We can 
do this most clearly by making a side-by-side comparison of the regularization schemes in 
Table | using the error analysis advocated by Lepage ||. This comparison clearly illustrates 
which schemes behave like a true EFT. Those familiar with the successful phenomenological 
models for NN scattering such as the Reid ||15|| , Bonn [Q, and Paris potentials, which 
already reproduce the data well, may wonder why there is a need to reformulate the problem 
in terms of an EFT. We use the Reid potential as an example to compare conventional and 
EFT approaches. 

We start by reviewing the expected behavior and radius of convergence of an EFT and 
outline the various regularization schemes in Section II. We then reiterate Lepage's discussion 
on how to analyze effective field theory behavior in terms of error plots || in Section III. We 
also show how to numerically extend his approach beyond second order in the momentum 
expansion and how to determine the radius of convergence graphically from the error plots. 
This procedure allows us to compare different regularization methods on the same footing. 

Our experience with a variety of different models suggests that the issues mentioned above 
are generic. Therefore we use the simple delta-shell potential without pion contributions as a 



TABLE I. Regularization schemes and their abbreviations used throughout this paper. 



name 



regularization scheme 



CR[G] 
DR[MS] 
DR[PDS] 
dibaryon 



Cutoff regularization with gaussian weighting in the potential || 
Dimensional regularization with modified minimal subtraction 



Dimensional regularization with power divergence subtraction 12] 

Additional low-energy degree of freedom associated with a bound or nearly bound 

state |JlC]. 
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convenient "laboratory" for comparing regularization methods in Section IV. We find that all 
of the regularization schemes considered in this paper have a proper radius of convergence 
except for dimensional regularization with modified minimal subtraction (DR[MSJ) when 
applied as in Ref . @ . Furthermore, a good radius of convergence is linked to the naturalness 
of the constants from the effective lagrangian and the ability to make reliable error estimates. 

An effective field theory should provide a low-energy realization of the S-matrix. This 
implies that after fitting the potential to a given order by the (scattering) phase shifts, 
estimates of the bound state energies should have similar error scaling. We present results 
for the bound state in Section IV and our conclusions in Section V. 



II. DEFINING AN EFFECTIVE FIELD THEORY 

The strong interaction at low energy can be described by an effective lagrangian of 
hadrons, whose form is constrained by the symmetries of the underlying theory, QCD. 
Restricting the characteristic momentum p of the interaction to be less than a scale A 
implies that physics at larger mass scales, such as from the exchange of heavier particles, 
is not resolved. This separates the physics into a long-distance part, given by the light 
dynamical fields in the effective lagrangian, and a short-distance part consisting of the 
heavy degrees of freedom and entering only through renormalization of effective lagrangian 
couplings. An exception is the nucleon, which due to baryon number conservation can be 
present as a heavy source. 

The pion is the lightest of the hadrons and a prime candidate to be included as a 
dynamical degree of freedom. Its Goldstone nature leads to a systematic treatment, known as 
chiral perturbation theory ||16|| , that has been thoroughly studied. No systematic treatment 



is known for non-Goldstone particles such as the p meson. Its exclusion from the effective 
lagrangian sets the scale of new, underlying physics A ~ m p . 

We stress that the effects of the short- distance physics on the long-distance (low- 
momentum) observables can be reproduced by generic terms in the effective lagrangian, 
organized in powers of derivatives over A, or a momentum expansion [[17]]. The coefficients 



of these terms summarize the remnants of the short-distance physics. In the future, these 
renormalized coefficients could in principle be calculated on the lattice. Until then, they 
must be fixed by existing data. 

As each additional coefficient in the lagrangian is fixed, the error in the S-matrix should 
improve by a power of p 2 /A 2 [[]]. (If the interaction contains single powers of momentum 
coupled, for example, to an external field, the improvement may be by powers of p/A.) 
This identification of the accuracy with powers of A, called power counting, implies that the 
calculation of every observable from the same EFT will give the same systematic momentum 
dependence on the error, as long as the corresponding effective operator for each observable 
is determined to the same order in A. Furthermore, the dimensionless coefficients of the 
terms in the effective lagrangian are typically of order unity (natural) when the cutoff is 
chosen on the order of the physics not contained in the lagrangian. All the above features 
are favorable for systematic predictions and an assessment of the corrections. 

At best, this expansion breaks down when the momenta in the processes involved are 
comparable to A, at which point the short- distance physics begins to be resolved. This point 
is referred to as the radius of convergence of the EFT. We will see below (e.g., the second 
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plot of Fig. |l]) that this also has a clear graphical interpretation as a convergence of the error 
involved when an increasing number of orders in p/A are taken into account in the effective 
field theory. 

Systems consisting of two (or more) nucleons with momenta well below the scale A can 
be treated nonrelativistically. Interactions between the nucleons lead to infrared divergences 
that require a summation to all orders of a certain class of diagrams This can be achieved 
by solving the Schrodinger equation or, equivalently, the Lippmann-Schwinger equation 

T = V + VG T . (1) 

The short- distance effective lagrangian is schematically mapped to a coordinate-space poten- 
tial consisting of highly singular terms with delta functions and derivatives of delta functions. 
Only after specifying a regularization scheme is the effective potential V defined. The energy 
dependence of V can be traded for three-momentum dependence by using the equations of 
motion JIB |. 



The different ways the potential can be regularized lead to different candidate effective 
field theories. In perturbative applications the predictions of an EFT are independent of the 
regularization method. However, this result has not been established for nonperturbative 
applications; indeed, some doubts have been raised about the equivalence [f^Pl. 

Cutoff regularization is a physically intuitive method for dealing with divergences in 
which momenta greater than the scale of new physics A are explicitly suppressed. This can 
be implemented by regulating the momentum integrals or by regulating the potential 
itself 0. We focus on the latter approach and simplify our discussion by concentrating on 
S'-wave scattering. We assume the long-distance physics is properly taken into account in the 
potential and focus only on the short-distance EFT potential. After defining a regularized 
delta-function: 5„(r) = e~ r l 2a /(27r) 3//2 a 3 , the 5-wave potential can be taken as 

(r'|F CR[G] |r) = Ana 2 5 3 (r - r') (c5 3 (r) + \da 2 V 2 ^(r) + \ea A V 4 5 3 (r) + ...), (2) 

which is equivalent to including a gaussian suppression of the transferred momentum q = 
p — p', as seen from the Fourier transform: 

(p'l^Gilp) = Vor[G](p, P') = 4 ™ 2 ( c - \ d <? a * + \ e ^ + • • •) eWaV2 ■ ( 3 ) 

We will refer to this method as Cutoff Regularization with Gaussian weighting, or CR[G] 
for short. The coefficients c, d, . . ., are dimensionless. The factor of 1/4tt picked up by 
each additional term in the Born series requires a Air to be factored out in order to render 
these coefficients natural [|J, i.e. of order unity, as we discuss further in Section IV. The 
coefficients are determined order-by-order from matching to the available data. To fit the 
first n constants in the potential requires at least n points of data. This data should be 
taken at as low momentum as feasible to minimize the contribution from terms of 0(q 2n a 2n ) 
that have been omitted. 

We need to make a few comments on the simplicity of Eqs. ([|) and (|3|). First, terms 
that do not contribute to S-w&ve states such as p • p' are included in the potential (as can 
be seen by expanding the q 2 term for example) to simplify the position-space expression 
that we use to evaluate the amplitude ||. However, the regularized delta-function mixes 
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these terms and they end up contributing to S-w&ve scattering at higher orders in the a 2 
expansion. The net effect to the order we are working is to shift the constant e by a natural 
amount. Since we are only interested in the naturalness of the constants in this paper and 
not their exact value, we will not discuss the effect of such contributions below. 

Second, terms such as d' Vf^(r) -V and e' a 4 (V 2 <^(r)V 2 + /i.c.) are omitted from Eq. (0). 
These terms, when applied with on-shell solutions of the Lippmann-Schwinger equation, only 
serve to modify the coefficients already taken into account in Eq. (|2|) JIDfl . Since we fit at 
most three constants below, we can take d! — e' — without any loss of generality to our 
discussion here. Extending the analysis off-shell and to channels other than S'-wave requires 
a more careful analysis of these contributions. 

One advantage to the cutoff method is that an increase of the momentum cutoff 1/a 
beyond the scale for new physics A is signaled by unnaturally large coefficients in Eq. ([3]) 
and degraded EFT estimates for the scattering amplitude. This behavior was illustrated by 
Lepage || and we have verified his results in our analysis. It provides a way to determine 
A for an unknown potential or data and suggests 1/a ~ A as a natural value for the cutoff. 
This choice ensures that all the physics treated correctly at low-momentum is still taken 
into account, but the higher- momentum physics is suppressed. Note that while the cutoff 
is roughly fixed by the physics, it is not fine-tuned to data as is done for the coefficients in 
Vcr[g] 0- Instead, the cutoff is a generic parameter of the EFT that dictates the radius of 
convergence and naturalness of the constants. 

Some disadvantages of the cutoff method are that the nonperturbative solution of the 
Lippmann-Schwinger equation requires numerical techniques and does not transparently 
exhibit simple power counting. Furthermore, both chiral and gauge symmetries are broken 
with a cutoff and require additional counterterms to restore them. 

An alternative approach is to use dimensional regularization, which preserves the symme- 
tries of the underlying theory and analytically can be shown to have simple power counting. 
This can be implemented for on-shell solutions of Eq. ([I]) by using the potential 



and using the modified minimal subtraction prescription DR[MS] or power divergence 
subtraction DR[PDS] [|1^] on the momentum integrals. Here A s is a scale introduced to make 
the coefficients dimensionless. Since DR[MS] has no scale (such as a cutoff) associated with 
the divergences, the natural size for A s after fitting to the data is dictated by physical low- 
energy scales such as the scattering length a s and effective range r e . As shown in Refs. [FU§, 
keeping the nonperturbative amplitude to all orders in p 2 , the momentum expansion is in 
powers of p 2 a s r e /2, which breaks down at very small momentum when the scattering length 
is large. This is a critical issue, since the AW system is known experimentally to have an 
almost bound state in the l So channel and a weakly bound state in the 3 Si channel (the 
deuteron). These are reflected in a large scattering length corresponding to momentum 
scales of about 35 MeV and 90 MeV respectively. 

One way to fix this pathology of DR[MS] is to introduce a low-energy degree of freedom, 
referred to as the dibaryon, into the lagrangian to parameterize the rapid energy dependence 
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of the amplitude from the bound or almost bound state |9|].[] Taking the large scattering 
length explicitly into account removes it from the residual momentum expansion and im- 
proves the radius of convergence. Another solution is to use DR[PDS], which allows an 
arbitrary scale \x to be included in the subtraction. This can be shown to give a momen- 
tum expansion in p 2 r e /2(fi — l/a s ), which is well behaved for large scattering length if an 
appropriate /i is chosen fT2| . With the power counting scheme advocated in fl~2|1 , DR[PDS] 
can be shown to produce a fi independent result. Since [i = corresponds to DR[MS], this 
means a specific power counting scheme cures the problems of dimensional regularization as 
noted in [|3,r0|. For comparison, all references to DR[MS] below refer to the naive solution 



of the Lippmann-Schwinger equation, and DR[PDS] will refer to this modified power count- 
ing scheme. All three dimensional regularization techniques can be solved analytically for 
short-range potentials. 

Finally, we contrast conventional nonrelativistic NN phenomenology to the effective 
field theory approach. The most important difference is that while the phenomenological 
potentials incorporate basic pion physics such as one-pion exchange, none of them has a 
systematic and complete inclusion of the long-distance Goldstone boson physics. But what 
about the systematic incorporation of short-range physics, as we consider here? If we take a 



Reid-style potential []IS| as an example we can see similarities. The original Reid potential 
consists of a sum of Yukawa interactions, with fixed masses chosen as integer multiples of 
the pion mass, and coefficients that can be varied to fit the data. In general for S-w&ves, 

^Reid(p, P') = 2 + ~~ 2~~T 2 + ■■■ , (5) 

we have a series of Yukawa potentials that can be viewed as regularized delta functions with 
different regulator masses mj. Since we deal only with the short range interactions, we are 
interested in setting these masses on the order of the scale of underlying physics instead of 
m v . If the terms in Eq. (||) are combined, they generate a Taylor expansion in q 2 multiplied 
by a function that suppresses large momentum, similar to the effective field theory potential 
Eq. (0). However, since each term has a different mass, a clear separation scale A is not 
identified. 

Furthermore, the conventional procedure for fitting the constants has very different con- 
sequences. Reid sought the best global fit, varying the constants to minimize the \ 2 fit to 
the data over the entire range of momentum considered. This means he used a weighting 
based only on the errors in the data. Although this optimizes the overall agreement with 
experiment, it neglects the fact emphasized in effective field theories that a theoretical error 
of order p 2n /\ 2n is also present in the short-distance physics. 

The systematics of an EFT are obscured or lost completely unless this theoretical er- 
ror is also taken into account in the weighting of the fit. The ability to relate observables 
and processes to each other depends on the ordering of corrections in a well-defined mo- 
mentum expansion. Furthermore, only in such a context would constants larger or smaller 



1 This can also be shown to arise from a certain resummation of terms in any regularization scheme 
[E.IC]. For DR[MS], however, it is a requirement for a proper radius of convergence as we will see 
in Section IV. 
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than expected be a signal of new physics or symmetry constraints when matching to QCD. 
Therefore, the predictability of the effective field theory is intimately connected with a con- 
sistent treatment of the error, which requires both the empirical and theoretical errors to be 
taken into account. This is not considered in many EFT analyses and so we summarize the 
philosophy of effective field theory in the next section and show how to organize the analysis 
to reap the full benefits. 

III. PHILOSOPHY OF EFT ANALYSIS 

In this section, we describe how to analyze and compare the different regularization 
schemes in a manner commensurate with the features of an effective field theory outlined in 
the previous section^]. The potential is not a measurable quantity, so instead we must use a 
scattering observable such as the phase shift 8(p) to determine the constants. Inserting the 
potential Eq. @ into the Lippmann-Schwinger equation Eq. ([I]) and specifying a regular- 
ization scheme, we can solve for the amplitude and determine the phase shift by one of the 
following two equivalent relations 

In the left plot of Fig. [[[ we show the results from using CR[G] to fit one, two, and three 
constants in Eq. (|3]), as compared to the exact S- wave phase-shift for the delta-shell potential 
PHI , which models the underlying physics and is discussed in detail in Section IV. A first 
glance at the plot shows the approximation to the phase shift improves as more constants 
in the potential are fit. However, a second look shows that it is very difficult to gather 
any quantitative information from the plot. At what point the curve deviates enough to be 
considered inaccurate is not clear, and the radius of convergence of the EFT expansion is 
completely obscure. 

It is therefore more informative to plot not the phase shift itself, but the error in the 
phase shift: |^ e fr — ^truel @- If the effective field theory follows proper power counting, then 
this error should improve by two powers of momentum as each additional coefficient in the 
potential is fixed [see Eq. (|4])]. However, a simple calculation shows that every term in the 
momentum expansion of the phase shift contains the scattering length JR. This means that 
the error in the phase shift could be numerically sensitive to a large scattering length and 
contaminate the power counting. 

This is not an inherent problem. We avoid this issue with no loss of generality by plotting 
the error of p cot 8{p) instead: 

A[pcot5(p)] = \p cot 8 cS (p) -pcot5 truc (p)| . (7) 

Since a large scattering length is synonymous with a near bound state (or pole) in the 
amplitude, Eq. (@) shows that this pole is cleanly mapped only to the first constant in a 
momentum expansion of pcot8{p). It is known from conventional scattering theory that 



2 This analysis builds on the comprehensive discussion in Ref. ||. 
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FIG. 1. The phase shift 5{p) (left) and the error in pcot 6(p) (right), each plotted as a function 
of p for the delta-shell potential with a weakly bound state, as discussed in Section IV. The solid 
line is the exact result and the dashed lines show the CR[G] fit for one, two, and three constants. 

this combination has a well denned expansion in p 2 for short-range potentials known as the 
"effective range expansion" 

1 i 

p cot 6(p) = h -r e p 2 + v 2 p 4 + • • • , (8) 

CL S Z 

which defines the scattering length a s and effective range r e . When long-range potentials are 
included, Eq. ([8]) is only valid at low momentum or is even inapplicable. In this case, one 
must define a modified effective range expansion ||21|| . Since the effective theory contains the 
same long-distance physics as the true underlying theory, Eq. (0) can be modified to have 



a clean momentum expansion [fl9 |. For short range potentials as considered here, it suffices 
to use Eq. (0). 

This now gives a rigorous way to fit the coefficients in the effective potential V(p, p') 
to data. After evaluating the combination pcot5 e fr(p) in the effective theory, we subtract 
it from the true result (either data or an exact solution to a model problem) and fit the 
difference, 

2 4 

Apcot5(p) = a + P^ + ^ + ... , (9) 

to a polynomial in p 2 /A 2 to as high an order as possible. The convergence rate of Eq. (|9D 
determines the radius of convergence of the EFT, so we use this observable in our anal- 
ysis below. The use of other observables, such as the bound-state energy, is discussed in 
Section IV. 

Using a spread of momentum near zero, the polynomial fit should be weighted with both 
the expected theoretical error in momentum and any additional experimental noise. The 
resulting coefficients a, [3, 7, . . . are then minimized with respect to variations in the effective 
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potential constants c,d,e,... using an optimization code. In practice, this method is more 
robust and numerically stable than matching the values of p cot 5(p) at discrete points to fix 
the constants. This allows us to extend the analysis of Lepage || beyond second order as 
shown in the next section. We also note that such a procedure is needed when matching to 
data even when the EFT observables can be calculated analytically. 

The number of coefficients that can be minimized is given by the number of constants 



retained in the effective potential. We used DPOLFT from package SLATEC \T2\ to find 
the polynomial fit and MINF |23[| , which is based on the Powell method, to carry out the 
minimization. Normal accuracies in minimization using double precision numbers with this 
method are 10~ 15 or better. 

Plotting the error in p cot S (p) as a function of p on a log-log plot, we expect a straight line 
with slope given by the dominant (lowest) power of p/A in the error |6|]. As we include more 
constants, the slope in this error should increase, signifying the removal of higher powers of 
p/A. The second plot of Fig. |I] clearly demonstrates the order- by-order improvement in the 
amplitude as more constants are added to the effective potential. With one constant the 
slope of the error is two [i.e., 0(p 2 /A 2 )} and increases by two with each additional constant. 

This plot also gives a clear graphical interpretation of the radius of convergence of the 
EFT as the point where the lines of error converge. The second plot of Fig. [I] shows that 
A ~ 1 GeV. For this purpose, the error plot is much more informative than plots of the 
phase shift itself. Once the momentum is on the order of the cutoff, the effective theory 
breaks down as the short-distance physics is resolved. Only by inclusion of physics at and 
above the scale A can the effective field theory be taken beyond this point. We will discuss 
further consequences of this type of analysis in the next section. 



IV. ILLUSTRATION WITH THE DELTA-SHELL POTENTIAL 

In this section, we illustrate the points made above by using EFT techniques with the dif- 
ferent regularization methods to systematically describe the "unknown" short-range physics 
of a specific example. We could use NN scattering data, but for our purposes it is more 
convenient to use an exactly solvable potential to serve as data in order to have a clean 
understanding of what features are important. The delta-shell potential has been used in 
the past to simulate the large scattering length found in NN scattering [pD| |. Kaplan used 



this potential to illustrate the benefits of the dibaryon approach [|]. He found upon con- 
sidering NN scattering that the inclusion of long-distance pion physics did not change the 
conclusions. This agrees with the experience of other authors |j],[7j that the addition of pions 
as long-range interactions does not affect the properties of the short-range expansion. The 
delta-shell potential is therefore a sufficient model for our purposes here. 

The potential can be written in terms of the nucleon mass M, the coupling g, and the 
range of the potential r$. This short-range potential represents the new physics of our 
underlying theory, and so we take r = 1/A below, 

V^ e (r) = -g^d(r-l). (10) 



M V A 

It has exactly one bound state for g > 1 and no bound states for g < 1. Scattering with 
p > A probes the details of the potential, so we expect A to be the radius of convergence 
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of a well-tuned EFT. The scattering length becomes very large for g near 1, whereas the 
effective range (and the rest of the terms in the effective range expansion) are of natural size 
for all g: 

9 1 2(o + l)l 

•' ' (11) 



9-1 A' 3g A 

The scattering length in the 1 So channel of iViV-scattering can be modeled by choosing 
(g,A) = (0.99, m p ). This potential Eq. fllPf ) with different p's will be the "laboratory" from 
which we compare the different regularization schemes. Results for actual measured data 
and the inclusion of pions will be discussed in a future publication ||19|| . 

We will first gain some intuition by graphically reproducing Kaplan's result that DR[MS] 
has a small radius of convergence if the scattering length is large. From Eq. (|TTD , we see that 
choosing g = 0.99 gives a scattering length one hundred times larger than choosing g = —10. 
At the same time, we investigate the effect from the presence of a bound state by taking 
g = 1.01. The effective potential is given by Eq. (|J) with the mass scale A s associated with 
the inverse delta-shell radius and the prescription of using DR[MS] on all divergent integrals. 

As mentioned above, the momentum expansion for DR[MS] can be shown analytically to 
be p 2 a s r e /2 |7|,|9| . This implies using Eq. (|TT| ) that the radius of convergence for g = 1.01 and 
0.99 should be roughly 1/10 that of the g = — 10 case. Fixing the constants in the potential 
Eq. (|]) by matching p cot 5(p) as outlined in the previous section, we produce the results 
in Fig. |2|. The dashed lines show the DR[MS] results for one, two, and three constants 
respectively. Indeed all three lines converge to p/A ~ 1 for g = —10 and p/A ~ 0.1 for 
g = 1.01. The results for g = 0.99 fall on top of the g = 1.01 results and are therefore not 
shown. This implies the presence of a bound state as opposed to an almost bound state does 
not matter, but the size of the scattering length does. It would be difficult to draw these 
conclusions had we only plotted the phase shift itself (the left plot of Fig. p. Of course, we 
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TABLE II. Effective potential for g = — 10 (small scattering length) to three different orders 
for different regularization schemes. Aa = 1 for CR[G] and /i = A for DR[PDS]. 





c 


DR[MS] 

d e 


c 


DR[PDS] 

d e 


c 


CR[G] 
d 


e 


0{p 2 /A 2 ) 


0.758 




8.33 




1.49 






eV/A 4 ) 


0.758 


-0.206 


8.33 


-25.0 


2.67 


-1.04 




0(p 6 /A 6 ) 


0.758 


-0.206 0.0672 


8.33 


-25.0 38.2 


2.67 


-1.58 


3.75 



TABLE III. Effective potential for g = 1.01 (large scattering length) to three different orders 
for different regularization schemes. Aa = 1 for CR[G] and jx = A for DR[PDS]. 







DR[MS] 




DR[PDS] 




CR[G] 




c 


d e 


c 


d e 


c 


d e 


0(p 2 /A 2 ) 


84.2 




-0.842 




-1.42 




0(p 4 /A 4 ) 


84.2 


-5.64 x 10 3 


-0.842 


-0.564 


-0.946 


0.842 


0(//A 6 ) 


84.2 


-5.64 x 10 3 3.78 x 10 5 


-0.842 


-0.564 -0.152 


-0.937 


0.618 -0.181 



can also show these results analytically for this simple model, but the analysis applies much 
more generally, when part or all of the calculation is done numerically. 

The constants c, d, and e are given in Tables |I| and |ITT| . Their values are calculated 
numerically to at least 8 digits to produce the accuracy of Fig. |2] and agree with the analytical 
values. This serves as a check on our numerical routines. The correlation between the 
naturalness of the constants and the radius of convergence is apparent. The constants for 
the g = 1.01 case are extremely large, reflecting the breakdown of the effective field theory 
for DR[MS] much below the expected scale A. 

One way to fix this behavior is to introduce the dibaryon |8]-[T0|. This takes the large 
scattering length into account by explicitly introducing a low-energy s-channel degree of 
freedom into the effective lagrangian. For g > or g < — 1, the potential can be written as 

t/ ( <\ n y 2 n- 2n 2 _3ttA1 + (7 3A 2 l-g 

Wn(P,P)-C-^; C ~MA> y ~ ~NP- g ' A "2M^' {U) 



with E = p 2 /M always kept on-shell. Since there seem to be three constants fit in Eq. (|12|). 
one might think the dibaryon will have an error of 0(p 6 /A 6 ). However, the dibaryon am- 
plitude is only matched to second order in the momentum when deriving the relations in 
Eq. (|i~2l ) DH and indeed shows an error of 0(p A /A 4 ) for the two values of g in Fig. § (plot- 
ted as the dot-dashed line). The slope and magnitude of the error do not depend on the 
scattering length, as expected. 

We next repeat the calculations using the cutoff regularization method CR[G] with 1/a = 
A. There are various ways to numerically solve the Schrodinger equation with a cutoff, but 
we have found the following procedure to be particularly efficient and numerically robust. 
First, the variable phase method is used to solve for the phase shift. This is a differential 
equation, 
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FIG. 3. The error in p cot 5(p) plotted as a function of p/A for a small scattering length without 
a bound state g = — 10 and for a large scattering length with a bound state g = 1.01. 



S'( r ) = -— V(r) sin 2 Q»r + 5(r)) with 5(0) = 
V 



(13) 



which expresses the change in the phase shift as the potential is built up from zero at r = 
to its full value at r = oo. The boundary condition ensures that the full phase shift given by 
5(oo) is zero in the absence of a potential and defines the otherwise ambiguous multiple of it 
in the phase shift. We use the routine ODE from package ODE [22] to solve the differential 
equation. To obtain the accuracy shown in the plots, we needed to do a weighted polynomial 
fit of Ap cot S(p) up to p/A = 0.1. 

The solid lines in Fig. show that CR[G] does work regardless of the scattering length. 
In fact, with A ~ 1/a, the result is just as good as the dibaryon for the same number of 
constants. The values for the constants are given in Tables [TI] and [TTI|, showing they are all 
natural for both g's (although the third constant is somewhat small for g = 1.01). 

Note that as more constants are fixed, the lower-order constants are modified. This 
occurs because even after truncating the potential in Eq. (||) to a given order, it still contains 
all orders in p 2 from the gaussian factor. The nonperturbative solution of the Lippmann- 
Schwinger equation therefore can generate terms of any order in p 2 / A 2 . The amplitude itself 
is matched to the true result order-by-order in the momentum so the power counting of 
the potential is destroyed. This consequence of the cutoff regularization is not necessarily 
relevant since the potential is not an observable. It is interesting to note, however, that 
these modifications are relatively small. This is also true when adding the long distance 



physics of the pion in fitting to actual NN data [19 

In summary, Fig. |2] shows that all regularization schemes considered so far produce useful 
effective field theories for a small scattering length, but DR[MS] fails for large scattering 
length. A failure of the power counting in powers of p/A is reflected in unnatural constants 
in the potential. 

We now focus on the most recently proposed regularization scheme, DR[PDS] [|l^]. The 
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potential is given by Eq. (f|), but with a further subtraction of an arbitrary mass scale 
n/Air from the linearly divergent integrals. This term mimics the behavior of the cutoff 
method, although still allowing a simple enough form to be solved analytically. An additional 
prescription compared to the DR[MS] case is an expansion of observables to the same order 
in p 2 as the potential Eq. (f|). If this is not done, the results are /i dependent, with \i = 
reproducing DR[MS] and \x larger than the nucleon mass approaching the CR[G] result in 
Fig.0. 

Since we are only dealing with a short-range potential, the DR[PDS] prescription repro- 
duces the effective range expansion Eq. (|5p by construction. The DR[PDS] results in Fig. |3| 
are therefore \i independent, although the constants still depend on /i. We take \x = A to 
compare with CR[G]. This produces natural constants for g = 1.01 but somewhat unnatural 
ones for g = —10 as seen in Table [[]J. This is an accidental consequence of the momentum 
expansion being in powers of p 2 r e /2(/i — l/a s ), so that for g = —10 and /i = A [see Eq. (|ll~l)l 
the denominator is nearly zero. The best result in this case occurs for /i = 0, which repro- 
duces the DR[MS] result and has natural constants. This implies that the scale \x is not 
functionally equivalent to the cutoff 1/a in CR[G] since it does not always signal the onset 
of new physics at the scale A. However, the [i independence of the scattering length shows 
this prescription for power counting produces a satisfatory radius of convergence, even for 
DR[MS] (fi = 0) |Jig,|T| . 



For both large and small scattering length, DR[PDS] does quite well, with a radius of 
convergence p/A ~ 1. The CR[G] result is better for one constant since the cutoff generates 
an effective range r e close to the true result. Overall, DR[PDS] is a convenient method to 
produce reasonable analytical results, and depending on the problem at hand either CR[G] 
or DR[PDS] may be suitable. One should note, however, that only DR[PDS] provides a 
strict diagram by diagram power counting |T2| . 

Finally, we return to the Reid potential ||15|| . The original Reid analysis used a global 
fit and only one adjustable parameter in each channel, for a result with approximately the 
same error (roughly a few percent of the data) at all momenta. However, we can apply the 
EFT fitting procedure instead. If we use Yukawa masses comparable to A, we anticipate 
similar results to CR[G]. Indeed, if a low-momentum fit is done to the constant c±, the error 
plot is similar to the CR[G] result with one constant (Fig. Adding a second short-range 
Yukawa does as well as CR[G] with two constants, since the Yukawas play off each other to 
allow the next order error in q 2 to be removed. This interplay becomes increasingly complex 
at higher orders. Furthermore, the additional mass scales obscure (or smear out) the role of 
A as a scale that separates the known from the unknown physics in effective field theories. 
Traditional AW potentials such as Reid are well suited for global fits. Systematic predictions 
with controlled error estimates are more properly analyzed using an effective field theory. 

We now turn to an investigation of the binding energy. If the EFT is truly reproducing 
the S-matrix of the underlying theory order-by-order in a momentum expansion, it should 
reproduce the binding energies and other observables to the same order of accuracy as the 
phase shifts. We therefore use the binding energy prediction as a consistency check for our 
candidate effective field theories. 

We have already fit the potentials to a given order by the scattering phase shifts above, 
and we use these potentials without adjustment to solve for the binding energy. This can 
be done analytically for the DR schemes by finding the poles in the scattering amplitude. 
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FIG. 4. The error in the binding energy for the candidate effective field theories for representa- 
tive values of g from 1.01 to 2.5. The star signifies the absence of real binding energies for DR[MS] 
with two constants at the values of g > 1.25 considered. 



The exact binding energy -Ebind for the delta-shell potential is given by solving the equation 

1 = 1 - e- 2 ^ = VME^T d 

g 2r) ' ^ A 

There is only one bound state to predict in the delta-shell potential, and if it is shallow 
enough, even the effective range expansion with the values of a s and r e can determine its 
value. A better test is to increase g until E hin d is large and on the order of A, and use 
the EFT to determine the accuracy of the binding energy prediction as a function of this 
variation. If a true radius of convergence is present, the effective field theory should break 
down for .Ebmd/A ~ 1. The binding energy is 6.48 x lCT 2 MeV for g = 1.01 but quickly 
increases to 812 MeV for g = 2.5. 

We plot the relative error in the binding energy in Fig. |j. Both CR[G] and DR[PDS] 
show the clear power counting behavior and proper radius of convergence expected from a 
true effective field theory. This gives a graphical verification that the errors in the binding 
energy really do follow power counting rules. We have checked that the same behavior is 
seen when plotted as a function of the average momentum (p 2 ) / A. The dibaryon result 
also follows the expected error scaling. In contrast, the deficiencies of DR[MS] regularization 
seen for the phase shifts are manifested here as binding energies that do not follow the EFT 
error scaling, improving to a lesser degree than expected. In addition, for values of g > 1.25 
with two constants, the DR[MS] S-matrix shows no bound state with a real energy. 

Therefore, our results show that most regularization procedures demonstrate the charac- 
teristics of a systematic predictive effective field theory. The fit of more and more constants 
in the effective potential improves the predictive power order-by-order in the momentum 
expansion. The radius of convergence of the EFT is independent of the scattering length 
and is given by the scale where new physics enters. 
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V. CONCLUSIONS 



Attempts to apply effective field theory methods to the nonrelativistic riucleon-nucleon 
scattering problem and then to the nuclear many-body problem have been stalled because 
of controversies concerning the nature and limitations of an EFT expansion when used 
nonperturbatively. The familiar dictum that "calculated observables are independent of the 
regularization method" has been questioned in this context. To help resolve these issues we 
have made a direct comparison of the various regularization approaches. 

We have applied the error analysis suggested by Lepage to a model problem to compare 
cutoff regularization, two forms of dimensional regularization, and the dibaryon approach in 
the context of nonperturbative, nonrelativistic effective field theories. This analysis focuses 
on a key signature of EFT behavior: the systematic scaling of errors with momentum or 
energy. We summarize some points made by Lepage about applying cutoff effective field 
theory ||: 

• The cutoff potential is not an observable and is not amenable to power counting. 
Individual constants can change as higher orders are taken into account, but predictions 
for observables are still systematically improved. 

• Fits to data should be weighted by both the uncertainty in the data and the expected 
theoretical error from power counting. This applies to any regularization scheme. 

• The cutoff should not be taken to infinity but only roughly adjusted to minimize the 
error, which identifies the resolution scale of the underlying short-distance physics. 
Using the cutoff as a fine-tuned parameter is not as effective as the introduction of a 
new low-energy constant in the potential. 

Our results verify these points. We also observed that natural coefficients are correlated 
with an optimal radius of convergence. 

New numerical procedures set forth in this paper allow us to work to third order and 
beyond in the EFT expansion, which is necessary to obtain a clear graphical determination 
of the radius of convergence for a given observable. Such an analysis is required for a 
systematic fit to data regardless of the regularization scheme. 

We find that all of the regularization methods except for dimensional regularization 
with modified minimal subtraction are consistent with basic features expected from a useful 
effective field theory: 

• Each additional order in the potential leads to a systematic improvement in the am- 
plitude. 

• The radius of convergence for this improvement, when optimized, is dictated by the 
scale of new physics. 

• Other observables are predicted with the same accuracy as the amplitude at each 
successive improvement. 

Our results are consistent with the analysis of van Kolck that, with proper resummations, 
any effective field theory for short-range interactions is equivalent to an effective range 
expansion to the same order ||. 
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The CR[G] and DR[PDS] regularization schemes are each suitable for developing effec- 
tive field theories of many-nucleon systems. In future work we will use both schemes in 
extending our fitting procedure and error analysis to iViV-scattering (including pions and 
other channels) and then to nuclear matter. 
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